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Abstract 


Calculated dynamic blade loads are compared 
with measured loads over a range of yaw stiffnesses 
of the DOE/NASA Mod-0 wind turbine to evaluate the 
performance of two versions of the MOSTAS computer 
code. The first version uses a time-averaged coef- 
ficient approximation in conjunction with a multi- 
blade coordinate transformation for two-bladed 
rotors to solve the equations of motion by stan- 
dard eigenanalysis . The results obtained with this 
approximate analysis do not agree with dynamic 
blade load amplifications at or close to resonance 
conditions. The results of the second version, 
which accounts for periodic coefficients while 
solving the equations by a time history integration, 
compare well with the measured data. To explain 
the deficiencies of the first version and to exam- 
ine the validity of the transformation, an investi- 
gation was made with the aid of a hypothetical 
three-degree of freedom dynamic model. The exact 
equations of motion of this model were solved using 
the Floquet-Lipunov method both before and after 
applying the transformation. Next, the equations 
with time-averaged coefficients after applying the 
transformation were solved by standard eigenanaly- 
sis. It was found that the transformation is valid 
but the associated time-averaged coefficient ap- 
proximation is inadequate for dynamic analysis of 
two-bladed rotors. 
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y 


multiblade coordinates of hypothetical 
model 

power train damping constants 
power train damping constants 

blade mass moment of inertia of hypo- 
thetical model 

mass moment of inertias of pod 
power train torsional inertias 
power train stiffnesses 

integer, 1, 2, 3, ... 
equivalent linear spring of tower in 
the lateral direction of hypotheti- 
cal model 

equivalent rotational spring of blade 
in the edgewise direction 
blade mass of hypothetical model 
hypothetical model effective mass mov- 
ing in the lateral direction which 
includes tower, bed plate, and gen- 
erator 

M + 2M 
e 

flatwise bending moment 
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m 

blade mass per unit length of hypo- 
thetical model 

n 

integer, 1, 2, 3, ... 

q 

tower lateral degree of freedom of 
hypothetical model 

R 

radius of hypothetical model 

r 

hypothetical model blade coordinate 

S 

blade mass moment of hypothetical 
model 

s 

Laplace variable 

X,Y 

coordinate system of hypothetical 
model 

Xp.Xp.Xp 

coordinate system of pod 

r g 

generator torque 

Y r 

torque applied to the rotor 

6M 

y 

flatwise cyclic bending moment 

r 

’’o 

multiblade coordinate of hypothetical 
model 

5 1 ,5 2 

lead-lag degrees of freedom of blades 
1 and 2 of hypothetical model 

a 

standard deviation 

♦ lf * 2 ,...* 5 

power train degrees of freedom 

*n 

blade phase angle of nth harmonic 

<l< 

azimuth angle 

“l ,u 2’“3 

coupled frequencies of hypothetical 
model in the fixed coordinated sys- 
tem 

“f *“f 

frequencies of Floquet transition 

1 2 *3 

matrix 

u t 

tower lateral frequency of hypotheti- 
cal model 


blade lead-lag frequency of hypothet- 
ical model 

0 

rotational speed of rotor 

C) 

time derviative 

( )’ 

derivative with respect to ^ 

[ ] 

square matrix 

{ } 

column matrix 

[ ] T 

transpose of a matrix 


Introduction 


The analytical calculation of rotor blade 
loads is one of the most difficult problems in both 
rotary wing technology and horizontal axis wind 
turbine technology, since it involves solving a 
highly nonlinear dynamic response problem. In 
Ref. 1 several analytical methods for calculating 
loads on a hypothetical helicopter rotor were com- 
pared. In Ref. 2 seven computer codes for calcu- 
lating dynamic loads in horizontal axis wind tur- 
bines were compared on the basis of calculated 
loads, with steady state measured data as a stan- 
dard. One of the codes evaluated in Ref. 2 is 
MOSTAS ,3 a code developed for NASA by Paragon 
Pacific, Inc. (PPI) to calculate dynamic loads and 
stability in a complete horizontal axis wind tur- 
bine system. That evaluation of MOSTAS was of a 
preliminary nature. 
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More recently the scope of MOSTAS has been ex- 
tended to account for time varying coefficients in 
the governing equations of motion of the wind tur- 
bine system. For descriptive purposes, the earlier 
version of MOSTAS will be designated herein as 
MOSTAS-A and the extended version as MOSTAS-B. 

The objective of this study is to conduct an 
indepth evaluation of the MOSTAS-A and -B codes to 
determine their suitability for predicting dynamic 
loads and instabilities in two-bladed wind turbines. 
This evaluation was accomplished by comparing blade 
loads measured on the DOE/NASA Mod-0 wind turbine 
against the load predictions of the two codes. The 
Mod-0 wind turbine, shown in Fig. 1, is an experi- 
mental 100 kW horizontal axis machine which was de- 
signed and built as a part of the DOE/NASA wind 
energy program and is now in operation at the NASA 
Plum Brook Station near Sandusky, Ohio.^ 1 

For the purpose of comparison, blade loads are 
chosen because of their importance in the design of 
large rotor systems. Since the dynamic coupling 
between the rotor and the tower and resulting dy- 
namic blade loads are sensitive to the variations 
in yaw stiffness, the loads for different values of 
yaw stiffness are compared. With the help of these 
results the performance of both versions of the 
MOSTAS code for predicting dynamic loads was evalu- 
ated. This evaluation raised some questions re- 
garding the validity of two major assumptions used 
for two-bladed rotors in MOSTAS-A. These assump- 
tions for calculating system frequencies, dynamic 
loads, and dynamic stability for two-bladed rotors 
are that the results are not affected either by the 
use of a multiblade coordinate transformation^ or 
by the use of time-averaged coefficient approxima- 
tion (approximate method). To assess the validity 
of these assumptions an additional study was per- 
formed with the aid of a simple hypothetical mathe- 
matical model. This paper presents results of both 
the load comparison study and the hypothetical model 
study. 


Description of MOSTAS-A and -B Codes 

The MOSTAS code is a collection of programs, 
beginning with a rotor code called MOSTAB-HFW in 
which the flexibility of each blade can be repre- 
sented by one to four modes. The support of the 
rotor is assumed to be rigid. Details of the -HFW 
code are given in Ref. 6. MOSTAS-HFW prints the 
key results of the trim-search process and also 
generates two disk or tape data files to be pro- 
cessed by subcodes PROCES and ROLIM. PROCES per- 
forms a harmonic analysis of the blade loads and 
ROLIM forms the rotor model from the single blade 
linear math model produced by -HFW. During this 
process ROLIM uses a multiblade coordinate trans- 
formation for a two-bladed rotor described in Ref. 5 
and synthesizes a rigorous linear rotor model in 
periodic coefficients. 

Unlike the multiblade coordinate transformation 
for rotors with three or more blades, 7 the trans- 
formation for two-bladed rotors assumes a constraint 
relation between the multiblade coordinates. The 
effect of these transformations on the equations of 
motion is to remove from the coefficients the har- 
monic terms which are not the integer multiples of 
number of blades. 


The ROLIM linear rotor model is combined with 
linear models of other components of the wind tur- 
bine system to produce a coupled system model. The 
coupling program is called WINDLASS. WINDLASS has 
provision for five different elements of the wind 
turbine system: rotor, control system, power train, 

pod, and tower. In the process of coupling, all 
periodic terms in the coefficients in the rotor 
model are time-averaged over one period so that the 
resulting equations for stability can be solved by 
standard eigenanalysis . The corresponding method 
of solution is referred to herein as the approxi- 
mate method. 

The dynamic response of the coupled wind tur- 
bine system to harmonic loads is processed by two 
subcodes, DYNAM2 and REC0V2 . DYNAM2 calculates the 
frequency response to harmonic loads and REC0V2 cal- 
culates time history response. The total dynamic 
loads are obtained by adding MOSTAB-HFW fixed-shaft 
loads and perturbation loads from REC0V2. The ver- 
sion of the MOSTAS code obtained by combining the 
programs MOSTAB-HFW, ROLIM, WINDLASS, DYNAM2 , and 
REC0V2 is referred to as MOSTAS-A. 

The MOSTAS-A code was extended to include all 
periodic coefficients in the equations of motion and 
to handle the arbitrary time varying loads by adding 
another program called W1NDGUST ,® in which the equa- 
tions with periodic coefficients are integrated 
numerically. WINDGUST replaces the programs DYNAM2 
and REC0V2. From the time history response, the 
total dynamic loads are obtained by adding MOSTAS- 
HFW fixed-shaft loads and the perturbation loads. 
This version of the system code is referred to as 
MOSTAS-B. 


Description of Mod-0 Wind Turbine 

The Mod-0 wind turbine is shown in Fig. 1. The 
system consists of five main components which are 
the rotor, the control system, the power train, the 
pod, and the tower. The rotor operates at 40 rpm, 
generating 100 kW of electric power in an 18 mph 
wind. The two aluminum blades are attached through 
pitch bearings to a rigid hub which is connected by 
means of a low speed shaft to a 45:1 step-up gear 
box. A high speed shaft (1800 rpm) connects the 
gear box to the 100 kW synchronous generator. The 
power train assembly is contained within a pod which 
is mounted on top of a 93 foot steel truss tower. 

Yaw control allows alignment of the rotor axis with 
the wind direction. The rotor is located downwind 
of the tower. A brief description of the mathe- 
matical models of the various components used in the 
MOSTAS code is given in Appendix A. 


Results and Discussion 
Experimental Data 

Bending moments measured in the shank area of 
the Mod-0 blades (5% of span) are used to evaluate 
the load prediction-capabilities of the MOSTAS-A 
and MOSTAS-B codes. Selected load data are divided 
into two groups, as follows: 

1. Data Case I . Data Case I was defined in 
Ref. 2 for purposes of comparing and validating wind 
turbine computer codes. This case represents the 
response of the Mod-0 wind turbine in its initial 
configuration, with stairs in the tower and a single 
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yaw drive unit connecting the pod to the tower. 

The nominal wind speed at rotor hub height is 28 
mph. Complete time histories of load are available 
for this case. Data Case I is characterized by 
high blade loads and substantial rotor/tower inter- 
action. 


2. Variable yaw-stiffness data . The effect of 
yaw stiffness on blade loads was measured by chang- 
ing the yaw drive system to the following configu- 
rations: (1) free yaw, in which the yaw drive was 

disconnected completely; (2) single yaw drive, 

(3) dual yaw drive, consisting of two parallel sin- 
gle yaw drive systems; (4) dual yaw drive with pre- 
load, in which backlash and nonlinearity were re- 
moved from the dual yaw drive system by loading the 
two systems against one another to a level of about 
30,000 lb-in; and (5) fixed yaw drive, in which the 
pod was clamped to the tower by brakes. In all 
cases the wind speed was approximately 25 mph and 
there were no stairs in the tower, the latter hav- 
ing been removed permanently to reduce tower shadow 
loads on the blades. 


Blade Load Comparison for Data Case I 

Because of its pronounced response, Data 
Case I is chosen for an indepth comparison of ex- 
perimental and analytical data. Figure 2 shows 
time histories of experimental and analytical blade 
bending moments in both flatwise and edgewise di- 
rections for this case. The abscissa is the blade 
azimuth angle which is 0° and 360° when the blade 
points downward. The flatwise bending moment is 
positive if the blade bends into the wind; the 
edgewise bending moment is positive if the leading 
edge is in tension. The experimental time histo- 
ries are presented as a band whose upper and lower 
bounds enclose data for three consecutive cycles. 
These data were also presented in Ref. 2. 

The comparison in Fig. 2(a) shows that the 
maximum and minimum flatwise bending moments pre- 
dicted by both M0STAS-A and -B codes are almost 
equal. However, MOSTAS-A results show four peaks 
per cycle and those of MOSTAS-B show three peaks. 
Measured maximum and minimum values are smaller in 
magnitude than those of either code. Measured data 
show three peaks and the amplitude of these peaks 
decreases at a much faster rate than that predicted 
by either code. This conservative nature of the 
calculated flatwise loads may be the result of the 
following factors: (1) the assumption that the 

structural damping of the yaw drive system is zero, 
(2) the differences between the actual and analyti- 
cal values of tower shadow and wind shear, and (3) 
the nonlinear behavior of the actual yaw drive sys- 
tem. The first factor is probably the most signif- 
icant. Some efforts to account for the damping of 
the yaw drive system are being made. 

The edgewise load comparison in Fig. 2(b) 
shows that the time history of moment predicted by 
MOSTAS-A significantly differs from that of 
MOSTAS-B. The measured time history of moment is 
in agreement with that predicted by MOSTAS-B except 
for a slight difference in phase angle. The maxi- 
mum and minimum moments predicted by means of 
MOSTAS-A are smaller than those obtained with 
MOSTAS-B. 

Table 1 presents a comparison between measured 
and calculated data on the basis of maximum, mini- 


mum, steady and cyclic bending moments. These re- 
sults quantify the comparisons in Fig. 2. The 
steady (My) and cyclic (£My) flatwise bending mo- 
ments are defined by the following equations: 
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where My >max and My )m i n are the extreme values 
of flatwise bending moments during one revolution. 
Similar equations apply for edgewise bending mo- 
ments . 


Examination of Table 1 shows that correlation 
between measured and calculated bending moments is 
highly variable. Good correlation was achieved be- 
tween measured edgewise moments and MOSTAS-B pre- 
dictions. All other comparisons show generally 
poor correlation. 

To compare the measured and calculated data on 
the basis of harmonic content, the bending moments 
shown in Fig. 2 are expressed in the following 
Fourier Series: 

£M = V'' C sin(n* + t ) n = 0,1,2, .. . (3) 

/ . n n 

n 


where 6M is the cyclic moment load in lb-ft, C n 
is the amplitude of the nth harmonic in the lb-ft, 
is the azimuth angle of the blade, and v n is 
the phase angle of the nth harmonic. Table 2 con- 
tains measured and calculated values of harmonic 
data for Data Case I. The measured amplitudes and 
phase angles are taken from Ref. 2 and these are 
averages of the bounding values. For further com- 
parison purposes, each harmonic amplitude was nor- 
malized with respect to its cyclic load. These 
normalized measured and calculated values are com- 
pared in Fig. 3. The comparison shows that the 
measured harmonic amplitudes continually decrease 
with increasing number except the fourth and fifth 
harmonics in the edgewise direction and the third 
harmonic in the flatwise direction. The correla- 
tion between the measured and calculated harmonic 
amplitudes varies from very poor to very good de- 
pending upon the harmonic number. 

In summary, Data Case I load comparisons were 
inconclusive as to the validity of MOSTAS-A and -B 
codes. Correlation was highly variable, with 
MOSTAS-B correlation being some what better than 
that of MOSTAS-A. 


Blade Load Correlation for Variable Yaw Stiffness 


Because of the importance of the cyclic bending 
moments for predicting life of blades and because 
of the sensitivities of the blade loads to the var- 
iations in yaw stiffness, the experimental and ana- 
lytical cyclic bending moments are compared for 
different values of yaw stiffness in Fig. 4. The 
abscissa is the effective yaw stiffness which is 
the resultant of the torsional stiffnesses of the 
tower and the yaw drive system in series. The mea- 
sured median cyclic bending moments together with 
16th and 84th percentiles are shown in Fig. 4. The 
band between the 16th and 84th percentiles approxi- 
mates a tier band and contains the loads for 68 per- 
cent of the rotor revolutions . 
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The comparison of cyclic flatwise bending mo- 
ments in Fig. 4(a) shows that both the MOSTAS-A and 
-B calculated loads are generally larger than the 
measured loads. MOSTAS-B predicts dynamic load 
magnifications whenever a nonrotating system fre- 
quency is close to an even integer multiple of the 
rotor rotational frequency. These magnifications 
in loads should be present because the entire wind 
turbine system as observed in the nonrotating coor- 
dinate system is subjected to zeroth, second, 
fourth, and sixth harmonic loads generated by the 
rotor. In particular, the magnification predicted 
by MOSTAS-B due to 2/rev resonance involving pre- 
dominently yaw motion is in good agreement with the 
magnification observed in the Mod-0 wind turbine 
with a single yaw drive and a dual yaw drive with- 
out preload. Although the flatwise loads predicted 
by MOSTAS-A for single yaw and dual yaw drive 
stiffnesses are higher than those of experiment, 
the load variation with yaw stiffness close to 
these cases is not in agreement with measured 
values. Also, the MOSTAS-A code is unable to pre- 
dict expected dynamic load amplifications when the 
nonrotating system frequencies are 4/rev and 6/rev. 
The MOSTAS-B code does predict these resonances. 

The comparison of edgewise bending moments in 
Fig. 4(b) shows that both MOSTAS-A and -B code 
values are in agreement with the measured values 
for the free yaw, preloaded dual yaw, and fixed yaw 
cases in which the system is free from resonance. 
The edgewise moments predicted using the MOSTAS-A 
code for the single yaw drive and the dual yaw 
drive without preload are not in agreement with 
measured values. MOSTAS-B loads, however, are in 
good agreement with measured loads for these yaw 
stiffnesses near the 2/rev resonance. The lack of 
correlation between measured and MOSTAS-A cyclic 
flatwise and edgewise bending moments at or close 
to the 2/rev resonance clearly indicates deficien- 
cies in this version of the code. 


Evaluation of Codes for Dynamic Load Predictions 

Based on the comparisons in the preceeding 
sections, it is inferred that the MOSTAS-A code, 
which uses an approximate method (time-averaged 
coefficients) in conjunction with a multiblade 
coordinate transformation is inadequate for pre- 
dicting dynamic load amplifications in two-bladed 
rotors. The MOSTAS-B code, which also uses the 
multiblade coordinate transformation but accounts 
for periodic coefficients, predicts dynamic ampli- 
fications in agreement with measured values. Com- 
parison of the two codes and measured data on the 
basis of harmonic amplitudes is inconclusive. 

The apparent cause of the poor performance of 
the MOSTAS-A code is the time-averaged coefficient 
approximation, which was introduced so that eigen- 
analysis could be used to solve the system equa- 
tions of motion. The multiblade coordinate trans- 
formation, on the other hand, appears to be valid 
since it is also used in the MOSTAS-B code which 
was found to yield valid load predictions. In 
order to test these two hypotheses — that the time- 
averaged coefficient approximation is invalid for 
two-bladed rotors and that the multiblade trans- 
formation is valid — a study was conducted by means 
of a hypothetical two-bladed wind turbine model 
with three degrees of freedom. Details of this 
study are given in Appendix B. 


As shown in Appendix B, calculation of system 
frequencies by means of accepted Floquet 
theory? > 9 ,10,13 unaffected by the multiblade 
transformation. 'However, the time-averaged coeffi- 
cient approximation which is used in the MOSTAS-A 
code alters system frequencies so that they do not 
agree with the exact Floquet frequencies. In par- 
ticular, a condition of "ground resonance" in which 
two system frequencies are equal was predicted by 
the Floquet theory before and after applying the 
multiblade coordinate transformation but not by the 
approximate method. Thus, the hypothetical model 
study supports both hypotheses: (1) the multiblade 

coordinate transformation for two-bladed rotors is 
valid, and (2) the time-averaged coefficient approx- 
imation which is used in the MOSIAS-A code is in- 
valid . 


Conclusions 

An investigation was performed to evaluate the 
validity of the MOSTAS-A and -B computer codes for 
predicting dynamic loads and instabilities in two- 
bladed wind turbines. This investigation was con- 
ducted with the aid of Mod-0 wind turbine data and 
the theoretical dynamic behavior of a hypothetical 
wind turbine model with three degrees of freedom. 
The following principal conclusions are drawn from 
the results obtained: 

1. The MOSTAS-A code is invalid for predicting 
either loads or instabilities in two-bladed wind 
turbines. 

2. The MOSTAS-B code is valid for predicting 
dynamic loads in two-bladed wind turbines. Good 
correlation between calculated and measured blade 
bending loads can be expected for systems like the 
Mod-0 wind turbine generator. 

3. The cause of deficiencies in the MOSTAS-A 
code is the time-averaged coefficient approximation 
introduced to permit solution of the system equa- 
tions of motion by eigenanalysis . 

4. The multiblade transformation used in both 
codes and its constraint relation do not affect 
system frequencies and can be used for analyzing 
two-bladed rotors. 


Appendix A 

Component Models of Mod-0 Wind Turbine 

For the purpose of analysis, the Mod-0 wind 
turbine is divided into five main components. The 
mathematical models for these individual components 
are described below. 

1. Rotor model . The rotor has two metal blades, 
each 62.5 feet long and weighing 2000 pounds. The 
flexibility of each rotor blade is represented by a 
modal model which is defined as a series of blade 
mode shapes and frequencies along with a definition 
of mass properties. These properties are listed in 
Table 3. The blade aerodynamic forces are gener- 
ated using "quasi-steady" blade element theory. 

The effects of gravity, the wind shear caused by 
the boundary layer of the Earth, and the tower 
shadow caused by rotor/tower aerodynamic inter- 
ference are included. The wind shear is assumed 
to cause a 30 percent linear increase in wind speed 
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from the lowest to the highest point in the rotor 
disk. The tower shadow effect is considered by 
retarding the free-stream velocity in 15° sectors 
on either side of the tower center line. The 
amount of retardation used is 50 percent of nominal 
wind speed for the tower with stairs and 27 percent 
for tower without stairs, based on wind tunnel test 
data in Ref. 11. 

2. Control system model . The pitch control 
system model used in this study is of first order 
as shown in Fig. 5(a). Also included in the con- 
trol system is the model for the synchronous gen- 
erator, Fig. 5(b), which is represented by a spring 
and a damper. Through this model, the generator 
torque is applied to the power train and the reac- 
tion is applied to the pod. 

3. Power train model . The power train model is 
composed of a series of one-degree of freedom 
models, each representing a flexible shaft with 
rigid pinion and main gears at its ends. A sche- 
matic of the power train model of Mod-0 is shown in 
Fig. 6 and the properties of various components are 
listed in Table 4. 

4. Pod model . The pod model is essentially an 
interfacing device for the other wind turbine com- 
ponents. It is assumed to be an elastic body with 
a concentrated mass. The elastic pod model is made 
up of three beam elements, as shown in Fig. 7, and 
its mass is concentrated at grid point 2. The mass 
properties are given in Table 5(a). The symmetric 
stiffness matrix Kp of this model (of order 
24x24) was generated using standard engineering 
beam elements. The nonzero matrix elements on the 
diagonal and in the upper triangle are listed in 
Table 5(b). This stiffness matrix is compatible 
with an ordered vector (24x1) of pod grid point de- 
flections resolved to pod inertial axes, shown in 
Fig. 7. The first six elements of the deflection 
vector represent three translations and three rota- 
tions for grid point 1, the next six for grid point 
2, and so on. The pod also includes the yaw drive 
system which is represented by concentrated springs. 
The torsional stiffness of the yaw drive system is 
the major structural parameter which was varied 
during this study. The stiffness values used in 
the analytical calculations are listed in 

Table 5(c). The Mod-0 wind turbine yaw drive 
stiffness for the data cases addressed herein are 
given in Table 5(d). 


Appendix B 

Hypothetical Three-Degree Freedom Model 

The purpose of the hypothetical model, as 
stated earlier, is to assess the validity of the 
multiblade coordinate transformation for two-bladed 
rotors and the associated time-averaged coefficient 
approximation (approximate method) which is used in 
MOSTAS-A code. To this end, a very simple model 
with two blades is considered, in order, to mini- 
mize the complexity associated with the Floquet- 
Liapunov theory which will be used to solve the 
equations of motion with periodic coefficients. 

The aerodynamic and gravitional forces are not in- 
cluded. The degrees of freedom of this model are 
the tower lateral motion and blade lead-lag motion. 
The tower is represented by an "equivalent" mass and 
a linear spring, and each blade is assumed to be 
rigid. The flexibility of the blade in the lead-lag 
direction is concentrated at the hinge which is as- 
sumed to be on the axis of rotation. A schematic of 
this model is shown in Fig. 8. 


The development of the equations of motion for 
this model is a straight forward application of 
Hamilton's principle. 12 Without the derivation de- 
tails, these equations are as follows: 


M t iJ + K t q + S (; cos t - i 2 cos V)" “ 0 
1;^ + + Sq cos ^ = 0 > 

It 2 + ” S '^ cos ' v = 0 


(Bl) 


where 


M = /” R mdr 
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S *= J " R mrdr 




r R 2 
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M = M + 2M 
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5. Tower . The tower is represented by a series 
of mode shapes and frequencies along with a defini- 
tion of its mass properties. The model properties 
of the tower were generated using the finite ele- 
ment structural analysis program NASTRAN. In cal- 
culating the modes and frequencies of the tower as 
an isolated component, a concentrated mass and mass 
moment of inertia were placed at the top of the 
tower to approximate the mass properties of the 
pod, rotor, power train, and generator system. 

Then, when all the modules of the system were cou- 
pled together, the concentrated mass and mass mo- 
ment of inertia were removed. By this approach a 
better representation of the tower mass and elastic 
characteristics is achieved with fewer number of 
tower modes. A total of three modes, one lateral 
bending, one fore and aft bending, and one torsion 
mode, were used for the towe\ model. The proper- 
ties of these modes are liste in Table 6. 


and M e is the effective mass in the lateral di- 
rection which includes both the tower and the pod. 


For convenience, the above equations are non- 
dimensionalized and expressed in the state vector 
form 


where 


[P](u') - [ Q ] { u} - 0 


(B3) 


[P] - 


co* I 0 


S 

MR ’ 


0 
1 

0 0 


-s 
M t R 

0 0 

0 0 

1 0 

0 1 
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IQ) 


2 

0 

0 


<■> ■ {l 


0 

0 

0 

2 


‘i} 1 




(B7) must be equal when all the frequencies are ex- 
pressed in either the cixed or rotating-coordinate 
system. 

To check this equality in frequencies, 

Eqs. (B3) and (B7) were solved using the Floquet- 
Liapunov theory for various values of rotor speed. 
Since the frequencies obtained from this theory are 
indeterminant, one has to add an integer multiple 
of the basic frequency ' 7 to each frequency obtained 
from the Floquet Transition Matrix. ? Thus, the 
frequencies of Eq. (B3) and the frequencies of the 
Floquet Transition Matrix of Eq. (B3) are related 
(Bi) by the equation 


The multiblade coordinate transformation for 
two-bladed rotors was developed in Ref. 5 and it is 
given by the relation 


1 0 cos \ Sla 

0 1 -sin » cos f 

1 0 -cos > -sis • 

0 1 sis i -cos » 


(B5) 


vv 


The effect of this transformation on Eq. (B3) is 
the absence of first harmonic terms in the coeffi- 
cients and this effect is similar to the one exhib- 
ited by the multiblade transformation for rotors 
with three or more blades. 10 The above transforma- 
tion, however, assumes a constraint relation be- 
tween the multiblade coordinates a. and b and 
it is (from Eq. (B5)) 


a| cos 4 + cos 4 = 0 


(B 6 ) 


Combining the Eqs. (B3), (B4), and (B5), the re- 
sulting equations in the multiblade coordinates can 
be obtained and these are 

[A]{ ; ' } - [ B ] { ?} = 0 (B7) 


where 



k 


i * 1,2,3 
k = 1,2,3,... 


(B9) 


Where u 2> an< * ‘*■3 are th e r °tor collective, 

the rotor cyclic, and the tower frequencies, respec- 
tively; inf. are the corresponding frequencies of 
the Floquet Transition Matrix. The basic frequency 
for Eq. (B3) is 1/rev. The rotor frequencies in 
Eq. (B9) are in the rotating-coordinate system. 

For a later comparison, these rotor frequencies are 
expressed in the fixed-coordinate system and are 
tabulated together with the tower frequency in 
Table 7. It should be mentioned the system fre- 
quencies of Eq. (B3) are in agreement with the re- 
sults presented in Ref. 13. 


The frequencies of Eq. (B7) and the frequen- 
cies of the Floquet Transition Matrix of Eq. (B7) 
are related by the equation 


u . 

l 


f! 



n 


+ 2 k 


i = 1,2,3 
k = 1,2,3,... 


(BIO) 


The basic frequency for Eq. (B7) is 2/rev. The fre- 
quencies uj, u) 2 , and 03 of Eq. (B7) are in the 
fixed-coordinate system. These are also included 
in Table 7. 


1*1 


U> - {uj » 2 c„ 

0 0 0 


0 0 5*» • 1 “ 2 ‘ ST 


1 0 
0 1 


IB] - 


0 0 0 

0 0 -Y: 


(1 ♦ co« 2*)| 

0 
0 

0 
1 


0 

- sin 2% 

0 

0 




ti -0 




^ (» - y\) 


(M) 

The tower motion in Eqs. (B3) and (B7) is in 
the fixed-coordinate system, and the blade motions 
are in the rotating-coordinate system in Eq. (B3) 
and in the fixed-coordinate system in Eq. (B7). If 
the multiblade coordinate transformation is valid, 
the system frequencies obtained from Eqs. (B3) and 


Comparison of the frequencies of Eqs. (B3) and 
(B7) in Table 7 shows that the multiblade coordi- 
nate transformation does not effect the frequencies. 
Hence, the transformation can be used for two- 
bladed rotors. 

To check the validity of the time-averaged 
coefficient approximation used in MOSTAS-A code, 

Eq. (B7) after approximating the periodic coeffi- 
cients in Eq. (B 8 ) was solved by a standard eigen- 
analysis method. The resulting frequencies are 
also tabulated in Table 7 under the heading "Approx- 
imate method." 

The frequencies tabulated in Table 7 are also 
plotted in Fig. 9 where the variation of nondimen- 
sional system frequencies as a function nondimen- 
sional rotor speed is shown with nondimensional 
tower frequency as a parameter. Several interestr 
ing observations follow from the results. These 
are: ( 1 ) for rotational speed ratios be- 

tween 0.2 and 0.4, the agreement between the Floquet 
and approximate method frequencies is rather poor, 
( 2 ) for rotational speeds ratios between 2.2 and 
2.6 the Floquet theory predicts the possibility of a 
ground resonance instability but the approximate 
method does not. Ground resonance is a phenomena 
characterized by the coalescence of two system fre- 
quencies and by the positive value of the real part 
of the eigenvalue. The Mod-0 wind turbine operates 
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in the rotational speed ratio range in which the 
agreement between the Floquet and approximate method 
frequencies is poor. These observations clearly 
suggests that the approximate method which is used 
in MOSTAS-A is not adequate to predict either the 
coupled frequencies or the possibility of ground 
resonance phenomena and hence any other dynamic 
instabilities. 


13. Miller, R. , Dugundji, J., Chopra, I., Sheu, D., 
and Wendell, J., "Wind Energy Conversion, 
Volume III: Dynamics of Horizontal Axis Wind 

Turbines," MIT Aeroelastic and Structures Re- 
search Laboratory, ASRL TR-184-9. 
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TABLE 1. - COMPARISON OF MEASURED AND 
CALCULATED BLADE LOADS FOR DATA CASE I 


(a) Flatwise loads 


Component 

Bending moment, ft-lb 

Measured 

MOSTAS-A 

MOSTAS-B 

Max 

Min 

Steady 

Cyclic 

127 500 
-2 500 
62 500 
±65 000 

161 200 
-46 000 
57 600 
±103 600 

153 000 
-42 000 
55 500 
±97 500 


(b) Edgewise loads 



Max 

34 500 

19 700 

34 000 

Min 

-76 500 

-54 500 

-81 800 

Steady 

-21 000 

-17 400 

23 900 

Cyclic 

±55 500 

±37 100 

±57 900 



TABLE 2. - HARMONIC CONTENT OF DATA CAST I BENDING 


LOADS AT STATION 40 (5% SPAN) 


(a) Flatwise bending moment 



Experimental 

MOSTAS-A 

Amp, 

Phase, 

Amp , 

Phase, 

ft-lb 

deg 

ft-lb 

deg 

35 400 

0 

37 490 


31 200 

24 

35 920 


25 700 

23 

17 240 


17 400 

-38 

36 250 

20 

7 100 

-70 

9 848 

139 

7 800 

-113 

50 300 

147 

3 000 

-126 

1 649 

178 


MOSTAS-B 


Amp, Phase, 
ft-lb deg 


35 770 
38 030 
17 170 
51 150 
10 730 
15 470 
1 722 



(b) Edgewise bending moment 










































TABLE 3. - ROTOR BLADE PROPERTIES 







Mode 1 

Mode 2 

Spanwise 


Blade 

Distrib- 

Edgewise 

Flatwise 

Edgewise 

Flatwise 

location, 


twist, 

uted 

deflec- 

deflec- 

deflec- 

deflec- 

ft 

M 

rad 

mass, 

slug/ft 

tion, 

ft/rad 

tion, 

ft/rad 

tion, 

ft/rad 

tion, 

ft/rad 

0 

0 


0 

0 

0 

0 

0 

0 

3.33 

0 


0 

0 

.0024 

-.0111 

-.0415 

-.0064 

4.00 

4.5 


.556 

0 

.0045 

-.0210 

-.0777 

-.0120 

6.75 

4.5 


.556 

1.4593 

.0218 

-.0988 

-.3558 

-.0497 

15.63 

4.5 


.261 

.8627 

.1914 

-.7779 

-2.3437 

-.1966 

25.00 

3.925 

.119 

1.1071 

.5924 

-2.4824 

-6.3133 

-.5092 

30.83 

3.558 

.0715 

.9678 

.9457 

-4.1833 

-9.6614 

-.9097 

37.50 

3.138 

.0346 

.8087 

1.4202 

-6.7544 

-14.1638 

-1.6545 

46.88 

2.547 

0 

.6104 

2.1844 

-11.7445 

-21.5842 

-3.4899 

50.00 

2.35 

-.0087 

.5778 

2.4575 

-13.8582 

-24.3409 

-4.4038 

56.25 

1.956 

-.0233 

.3139 

3.0171 

-18.8374 

-30.2452 

-6.8061 

59.38 

1.759 


.2961 

3.2928 

-21.8434 

-33.3369 

-8.4313 

62.50 

1.563 


.2664 

3.5673 

-24.9446 

-36.457 

-10.1432 




frequency 

, rps 

11 

.98 

17. 

76 


TABLE A. - POWER TRAIN PROPERTIES 


Parameter 

Values 

Element no. 

J, slug-ft^ 

K, ft-lb/rad 

B, ft-lb sec/rad 

C, ft-lb sec/rad 

1 

615.8 

9.8*106 

388 

766 

2 

938. A 
9.8*106 
A80 
960 

3 

1.22 

5.55*10^ 

1.30 

2.6 

A 

1.A08 

2.36*10 4 

0.91 

1.82 























TABLE 5. - POD PROPERTIES FOR MOD-O WIND TURBINE 


(a) Pod mass and inertias 


Mass, 

slug 

I xx» „ 
slug-ft^ 

, I yyi ■> 

slug-ft 2 

*zz, , 
Slug-ft 2 

1X18 

5000 

78600 

78600 


(b) Nonzero values of symmetric pod stiffness matrix 


Value 


Matrix position - row. 

column 


2,8 

3,9 

8,14 

1 ! 9 

■ 1 


i WWSS. M 

2,2 

3,3 

14,14 




-6.0X108 

3,5 

3,11 

6,8 

9,17 


14,18 

6.0x108 

6,2 

5,9 

8,18 

11,15 

Bsm 


-7.5x108 

7,19 

8,20 





7.5x108 

19,19 

20,20 





9.9x108 

8,8 






-1.0x109 

1,7 

9,21 





1.0x109 

1,1 

13,13 

21,21 




1.24xl0 9 

9,9 






-1.5xl0 9 

8,10 

8,22 

11,19 

19,23 



1.5x109 

7,11 

7,23 

10,20 

20,22 



2.0x109 

5,11 

6,12 

10,22 

11,17 

11,23 

12,18 

2.75x109 

7,7 






4.0x109 

5,5 

6,6 

17,17 

18,18 

22,22 

23,23 

-l.OxlolO 

4,10 

7,13 

10,16 

12,24 



l.OxlO 10 

4,4 

16,16 

24,24 




1.2xl0 10 

11,11 






1.8xl0 10 

12,12 






2.4x1010 








(c) Yaw drive system stiffnesses 


Shear and compression 

1.0x108 lb/ ft 

Bending 

1.0x108 ft-lb/rad 

Torsional (Yaw) 

0-1.0x1012 ft-lb/rad 


(d) Yaw drive stiffnesses for experimental data cases 


Yaw drive configuration 

Free 

Single 

Dual 

(no preload) 

Dual 

(preload) 8 

Fixed 

Stiffness, ft-lb/rad 

0.0 

7.91x108 

15.81x106 

22.35x106 

l.oxio 12 


a 30 000 lb- in. 
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TABLE 7. - COMPARISON OF FREQUENCIES OF FLOQUET THEORY AND 
APPROXIMATE METHOD FOR HYPOTHETICAL MODEL 


Rotational 

speed 


Frequencies from Floquet theory 




Approximate method 
(Eq. (B7) with time- 

























Figure 1. - DOE/NASA 100 kW Mod-0 wind turbine (formally ERDA/NASA). 




BLADE AZIMUTH, deg 


(b) EDGEWISE. 

Figure 2. - Comparison of experimental and analytical 
bending moment time histories for single yaw case 
(at 5% blade span). 




AMPLITUDE 
CYCLIC LOAD 



(a) FLATWISE. 



HARMONIC NUMBER 
(b) EDGEWISE. 

Figure 3. - Comparison of measured and calculated 
harmonic contents of moment load cycles. 


CYCLIC EDGEWISE MOMENT. Ib-ft CYCLIC FLATWISE MOMENT. Ib-ft 


T 84 PERCENTILE ^ 

0 MEDIAN ^MEASURED 

1 16 PERCENTILE J 

MOSTAS-A 

MOSTAS-B 


YAW DRIVE CONFIGURATION 

a FREE 
b SINGLE 

c DUAL NO PRELOAD 
d DUAL PRELOAD TO 30 000 lb-in. 
e FIXED TO TOWER 


80xl0 3 2/ rev peak 




(b) EDGEWISE. 

Figure 4. - Comparison of experimental and 
analytical cyclic bending moment for the 
Mod-0 wind turbine (at 5% blade span). 


REFERENCE 

SPEED 


ROTOR 

SPEED 



BLADE PITCH ANGLf 


CONTROLLER ACTUATOR 

PROPORTIONAL GAIN, K p -0.0 
INTEGRAL GAIN, K r *0.03 sec 2 
TIME CONSTANT, t * 100.0 sec' 1 
LAPLACE VARIABLE, s 

(a) PITCH CONTROL SCHEMATIC. 


o 

<\i 


GENERATOR 

DISPLACEMENT 

GENERATOR 

SPEED 



GENERATOR 
TORQUE • 


EQUIVALENT STIFFNESS, K * 1260 ft-lb/rad 
EQUIVALENT DAMPING, C g - 20 ft-lb sec/rad 


lb) GENERATOR SCHEMATIC. 

Figure 5. - Control system model of Mod-0 wind turbine. 


ROTOR 


GEARBOX BELT GENERATOR 
DRIVE 


Tr 





Figure 6. - Schematic for power train model of Mod-0 wind turbine. 
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POD MASS 
INTERFACE 
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GEAR BOX AND 

GENERATOR 

INTERFACE 
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YAW DRIVE 
INTERFACE 



Y P 

ROTOR 

INTERFACE 


Figure 7. - Schematic for pod model of Mod-0 wind turbine. 



Figure 8. - Hypothetical three-degree of free- 
dom model (q is tower lateral motion. ^ 
and ^2 are lead-lag motions of blades 1 
and 2, respectively). 



FLOQUET THEORY 

APPROXIMATE METHOD 

1 * ROTOR COLLECTIVE MODE 

2 - ROTOR CYCLIC MODE 



ROTATIONAL SPEED, 

Figure 9. - Variation of frequencies with rotational speed of a 
three-degree of freedom hypothetical model - tower lateral 
to lead-lag frequency ratio, - 0.924, S/M^R ■ 0.0693, 
SR/I * 1.84. 
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